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We analyze the crossover from an /ic/e-periodicity of the persistent current in flux threaded 
clean metallic rings towards an /ic/2e-flux periodicity of the supercurrent upon entering the su- 
perconducting state. On the basis of a model calculation for a one-dimensional ring we identify 
the underlying mechanism, which balances the hc/e versus the hc/2e periodic components of the 
current density. When the ring circumference exceeds the coherence length of the superconductor, 
the flux dependence is strictly hc/2e periodic. Further, we develop a multi-channel model which 
reduces the Bogoliubov - de Gennes equations to a one-dimensional differential equation for the 
radial component of the wave function. The discretization of this differential equation introduces 
transverse channels, whose number scales with the thickness of the ring. The periodicity crossover 
is analyzed close the critical temperature. 
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I. INTRODUCTION 

Charged particles, which encircle a magnetic flux 
threaded region, acquire a geometric phase. This 
Aharonov-Bohm (AB) phase leads to quantum interfer- 
ence phenomena along multiply connected paths 1 . A par- 
ticular manifestation of the AB-effect is the persistent 
current in mesoscopic metal ringiP^, which is modulated 
periodically by the magnetic flux piercing the interior of 
the ring with the period of a flux quantum $o = hc/e for 
clean rings. 

Likewise, in superconducting rings the order parame- 
ter responds periodically to magnetic flux, as implied by 
the requirement of a single- valued superconducting wave 
function in the presence of a supercurrentP^SI. Measure- 
ments of magnetic flux trapped in small cylinders proved 
that the flux in superconductors is quantized in units of 



r7 ; 



The hc/2e superconducting flux quantum was 



corroborated by measurements of the hc/2e periodicity 
of the critical t empe rature of superconducting rings by 
Little and Parkg^^l, and by the hc/2e flux quantization 
of Abrikosov vortices^. 

The oscillations of the persistent current or the su- 
percurrent with respect to the magnetic flux implies 
the corresponding periodicity for all thermodynamic 
functions^. Two classes of condensate states have been 
identified which are not related by a gauge transforma- 
tion. In the thermodynamic limit, they are degenerate for 
integer and half integer flux values, which results in the 
observed < &o/2 periodicity. This degeneracy is however 
lifted for discrete systems^, which was implicitly under- 
stood already in the early works of Byers and Yang^ and 
by BrenigSS The lifting of the degeneracy can be made 
explicit through the evaluation of the supercurrent in suf- 
ficiently small rings^H Recently, n odal supe rconductors 
have been in the focus of research^ 15 * 16 * 17 * 18 -) as they al- 
low for striking differences in the excitation spectrum for 



flux sectors centered around integer and half integer $o 
values, respectively. 

For rings of s-wave superconductors, it is expected, 
that the hc/e periodicity is restored if. the rin g diameter 
is smaller than the coherence length^ 19 * 20 * 21 !. While the 
$0 and the $o/2 periods are well understood in metallic 
and superconducting rings, it has remained unaddressed 
how the periodicity evolves for such a small ring when 
the normal metal turns superconducting. Here we ana- 
lyze the periodicity crossover in a one-dimensional (1-D) 
model for a flux threaded ring at zero temperature, which 
allows for a transparent analytical treatment on the ba- 
sis of the Gor'kov equations in an external magnetic field 
(Sec.|n]). 

The gap equation for the current carrying supercon- 
ducting state is solved for finite size rings to evaluate the 
field dependence of the discrete energy spectrum and the 
supercurrent. We identify two components of the current 
with hc/e and hc/2e periodicity, respectively, whose mag- 
nitudes shift with the opening and increase of the energy 
gap in the superconducting state. When the coherence 
length of the superconducting ring is of the order of the 
ring size, only the ft.c/2e-periodic component remains. A 
similar analysis for the temperature driven crossover in 
clean and dirty ID rings has recently been published by 
Wei and GoldbariP. 

It is well known that a long-range ordered supercon- 
ducting state does not exist in ID. However, whereas 
thermal phase slips suppress a transition into the super- 
conducting state at finite temperature, quantum phase 
slips at zero temperature are rare events. Even if phase 
coherence is broken at certain instants in time and space, 
the supercurrent does not decay in the ring. We there- 
fore investigate only the zero temperature transition for 
the ID ring. 



Subsequently (Sec. Ill I we extend our analysis to rings 
of finite thickness ("annuli" ), which represent, from a for- 
mal point of view, multichannel systems. For the annuli 



2 



we present the periodicity crossover upon cooling through 
the superconducting transition temperature. For annular 
systems, that are confined to a 2D plane, we introduce a 
semi-analytical approach in which the numerical work is 
reduced to the solution of a one-dimensional differential 
equation for the radial component of the wave function. 
Its discretization allows to introduce a fixed number of 
(transverse) channels, whose number parameterizes the 
thickness of the annulus. 



II. ID RING 

We start from the tight binding form of the kinetic 
energy for a, ID ring with N sites as given by 



Ho 



(i) 



where the sum extends over all nearest-neighbor sites i 
and j; s =f, | denotes the spin, and t is the hopping 
matrix element. The vector potential A of an exter- 
nal magnetic field enters through the Peierls phase factor 
(fij = (e/hc) J? A • dr = 2iT(j)/N, where 4> = $/$o an d 
<j> is the magnetic flux through the ring. After Fourier 
transformation Ho becomes 



H - £ k~<k c l 



ks C ks 



(2) 



ks 



with the single-particle energy for a state with angular 
momentum hk 



ek-s = -2tcos 



k-(j) 
R 



(3) 



i? = N/2ir denotes the dimcnsionless radius of the ring 
and k = -N/2, N/2 - 1. If N is a multiple of 4, a fc- 
value exists for which e k = 0^, with two states exactly at 
the Fermi energy Ep — for = 0. To ensure a unique 
ground state, we choose \i = t/N which is placed in be- 
tween two single-particle energies e k - This is achieved 
for even N, which are not multiples of 4. All calculations 
were performed for this generic choice of N and fj,. 

The superconducting state in this strictly II? ring 
model is controlled by a BCS-typc Hamiltonian of the 
form 



H = Ho 



k . q 



[A* k (q)c^ 



where A k (q) is the superconducting order parameter for 
the formation of Cooper pairs with finite angular mo- 
mentum Hq and q G Z. The order parameter is ob- 
tained from the anomalous imaginary time Green's func- 
tion F{k,k',T- T ') = (T lCfcJ (r) C _ t1 (r')flby 



A k (q) = k B T^2J2 V ^' F ( k '' k ' 



(5) 
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FIG. 1: Energy dispersion of a ring with an order parame- 
ter A = 0.22t for 4, = (dashed line) and <j> = <f> c « 0.24t 
(solid line), where the indirect energy gap closes. The filled 
(empty) circles represent occupied (unoccupied) fc-states for 
a ring with N = 18. The asymmetry for ±fe scales with 



where w n = (2n — VjirksT is the fermionic Matsubara 
frequency for temperature T and Vkk 1 is the pairing in- 
teraction; T T is the time-ordering operator. 

A k (q) has to be determined self-consistently in the 
superconducting state. This is achieved by solving the 
equations of motion for the anomalous Green's func- 
tion and the single particle propagator G(k,r — r') = 
(T t ci~ s (t)c\ s (t')) , which is diagonal with respect to mo- 
mentum and spin. This leads to the self-consistent set of 
Gor'kov equations: 

G- 1 (k,Lo n ) = G ~ 1 (k,uj n ) 

+ Y,^k(q)G (-k + q,-cu n )A* k _ q ( q ), (6) 



F(k,k-q,uj n ) = G (k,uj n )A k (q)G(-k + q,-uj n ), (7) 

where Go(k, ui n ) — [ihuj n — e^.^]" 1 is the Green's function 
in the normal state. 

We assume that the unique ground state of the super- 
conducting condensate is characterized by a single integer 
quantum number q(<f). For rings larger than the coher- 
ence length, the g-number of the ground state advances 
to the next integer whenever crosses the flux values 
(2n - l)/4, n G Z, i.e., q{<j>) = floor (20 + 1/2), where 
floor(:c) is the largest integer smaller than x. Discrete- 
ness of the energy levels shifts the increment in q slightly 
according to the energy difference of even-q and odd-g 
states. Disregarding this small shift in a first approach 
(see comment at the end of this section), we take A k (q) 
of the form 



A k (x) =5(x-q(<f>))A k 



(8) 



For s-wave pairing, which is the only Cooper-pair state 
possible in a strictly ID system, A k = A is constant. 




= 26 and a self-consistently calculated order param- 
black lines: E-(k,(j>), green lines: E+(k,<f>). Upper 
"large gap"regime (V = 1.9t, Ax/ 2 ~ 0.30t); Lower 
'small gap"regime (V = Lit, A1/2 ~ 0.08t). Super- 
ivity occurs only in the odd-q sectors for V = Lit 
. [3| . The bold line marks the highest occupied state 
1. For the definition of A1/2 see text after Eq. (11 1. 



With this ansatz, we determine the Green's function from 
Eq. (j6j as 

G(k,uj n )-— TvTTT^ . B TZ W 



(ihui n — E + {k, <j))){ihuj n — E_(k, </))) ' 
where the two energy branches E±(k,q) are given by 



± VA 2 + e 2 (/c,0) (10) 



with e(k, 4>) = (e fc -*+e_A; +g _0)/2. The energies J5±(fc, 0) 
are plotted in Fig. IT as a function of k. The upper (E + ) 
and the lower branch are separated by an indirect 

energy gap, which closes at a critical value A c . For finite 
flux the dispersion is asymmetric with respect to an in- 
version in the angular momentum k — > ~k (see Fig. [TJ, 
and this asymmetry induces a finite supercurrent. In the 
"small gap" regime A < A c , both E + (k, <fi) and E-(k, 4>) 
can be positive or negative, whereas in the "large gap" 
regime A > A c , E + (k, 4>) > and E-(k, (j)) < for all k 
and cf> (see Fig. |2ja)). Close to Ep, E±(k,cf>) simplifies 
to 



E±(±k, 



q) ± VA 2 + l(t/R) 2 (11) 







FIG. 3: Solution of the self-consistency equation ( 13 1 for 
different values of the pairing energy V at T = 0. From top 
to bottom: V = 1.9t, 1M, 1.35t, Lit. 



where k > and I = 1 for even q and / = for odd 
q. The maximum direct energy gap in the even-g sectors 
is therefore An = a/A 2 + (t/R) 2 , whereas in the odd-g 



sectors it is A 1 = A. Eq. (11) shows that the shift of 
the eigenenergies scales with the ring size as 1/R in the 
"small gap" regime. 

By inserting G(k,uj n ) into the Gor'kov equation (J7|, 
one finds for the anomalous Green's function 



F(k, k - q,uj n ) 



(ihh) n - E + (k,4>))(ihio n - E_(k,4>)) 

(12) 

For a momentum independent pairing interaction Vkk' = 
V we obtain the self-consistency equation for A(0) from 
Eq.(|5]) by summation over uj n 



N 



E 



f(E-(k,4>))-f(E+(k,4>)) = j. 
2^A( ( />) 2 + e 2 (fc, ( />) V" 



(13) 



where f(E) denotes the Fermi distribution function. In- 
stead of lowering the temperature we explore below the 
transition into the superconducting state at zero temper- 
ature by increasing the pairing interaction strength V. 
The flux (j) affects the solution of the gap equation ( 13 1 



for A in two ways. For small-size rings the magnitude of 
A is mainly controlled by the energy of the level closest 
to Ep. If the quantity 6^ = min^ \e(k,<j>) — Ep\ > 0, a 



solution of Eq. ( 13 ) exists only above a threshold value 
of the pairing interaction. In the even-g sectors, this is 
the case for all values of <fi, whereas in all odd-g sectors a 
flux value <f> exists, for which 5^—0 and Eq. (13) has a 
solution for all V > (cf. Fig. a\. This is a consequence 
of the discreteness of the energy levels. In the strong 
coupling regime V 3> t, A is modulated only slightly by 
the flux. For weak coupling V ~ t, a solution Ai < A c 
is possible in the "small gap" regime, where Ai denotes 
the order parameter at half-integer flux values. In this 
case the energy gap closes at a critical flux <f> c in the odd- 
q sectors and E + (k, cj>) turns negative for the level closest 
to Ep (see Fig. [2]). Thus the dominant term in the sum 
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FIG. 4: Crossover from the /ic/e-periodic normal persistent current to the /ic/2e-periodic supercurrent in a ring with N = 26 
at T = 0. For this ring size A c w 0.24*. The discontinuities occur where the (^-derivative of the highest occupied state energy 
changes sign. 



of Eq. ( 13 1 switches sign and the solution for A vanishes 



discontinuously. This is equivalent to a breaking of the 
Cooper pair closest to Ep, which provides the main con- 
tribution to the condensation energy^MSJ. These features 
for the solution of the self-consistency equation are spe- 
cial for strictly ID rings. In these rings superconductivity 
is destroyed for velocities of circulating Cooper pairs ex- 
ceeding the Landau critical velocity, which is approached 
at (f, = </> c EH. 

With the discrete lattice gradient Vj/(i) = 
\ [f(i + 1) — f(i — 1)], the current is obtained from 



-te 



J{4>) = -rfyi-Vjffli-jy 



1=3 



(14) 



where n(k) — kgT^2 n G(k,uj n ) is the momentum dis- 
tribution function. The result is shown in Fig. [4j For 
V = A — one recovers the /ic/e-periodic saw-tooth pat- 
tern for the normal persistent current as discussed irP. 
With increasing A, new linear sections appear continu- 
ously. These are the sections where the order parameter 
is finite in the "small-gap" regime^. The occupied state 
closest to Ep contributes dominantly to the current, be- 
cause all other contributions tend to almost cancel in 
pairs. The discontinuities of the current occur where the 
(f> derivative of the energy of the highest occupied state 
switches sign (see Fig. [2]). These linear sections increase 
with increasing A; once they extend to a range hc/2e 
upon reaching the "large gap" regime, the current be- 
comes strictly ft,c/2e-periodic. 

We obtain further insight into the mechanisms, which 
determine the current periodicity, by analyzing A c . Ac- 
cording to Eq.([TT]), close to Ep, the maximum energy 
shift is t/(2R), and the condition for a direct energy gap 
(or E + (k, <fi) > for all k, <fi) and an /ic/2e-periodic cur- 
rent pattern is therefore A > A c = t/(2R). The corre- 
sponding critical ring radius is R c = t/(2A). 

It is instructive to compare R c with the BCS coherence 
length £ = hvp/(jrA), where Vp is the Fermi velocity 
and A the BCS order parameter at T = 0. On the lattice 



we identify vp — hkp/m with kp — ±7r/2a and m = 
h 2 /(2a 2 t); a is the lattice constant. Setting the length 
unit a = 1 we obtain £o = t/A and thus 2R C = £q. This 
signifies that the current response of a superconducting 
ring smaller than the coherence length is generally he/ e- 
periodic^l In these rings the Cooper pair wave function 
is delocalized around the ring. 

A second fundamental effect, which manifestly breaks 
the /ic/2e-periodicity, is the offset of the transition from 
even to odd center of mass angular momenta q with re- 
spect to evenly spaced flux values (2n — 1) hc/4e. This 
small offset was already observed in our previous numeri- 
cal evaluations for rf-wave loops 1 ' 1 . Vakaryuk 21 has traced 
this shift to the dependence of the internal energy of 
Cooper pairs on the center of mass state. For a BCS- 
model superconductor, this effect is fully incorporated in 
the Bogoliubov - de Gennes (BdG) evaluation of Ref. [TS] 
although the quasiparticle-like presentation introduces a 
different perspective. In the discussion in this section we 
disregarded the offset for the ID rings in order to focus 
on the aspects related to the opening of an indirect gap. 
In Sec. Ill we include the offset consistently in the BdG 
evaluation of the multi-channel annulus. 

It is worthwhile to note that the condition A > A c 
(or R > R c ) only refers to the periodicity of the super- 
current. It does not guarantee an /ic/2e-periodicity of 
the order parameter A or the total energy, but only of 
their derivatives. These quantities need a continuous en- 
ergy spectrum withdegeneracies for flux values which are 
multiples of hc/2^^. 



III. MULTICHANNEL RING: ANNULUS 

In this section we describe a superconducting loop of 
finite width as shown in Fig. [5] with an inner radius R\ 
and an outer radius i?2- For such an annulus, we choose 
a continuum approach on the basis of the Bogoliubov - 
de Gennes (BdG) equations. For integer and half-integer 
flux values, these equations can be solved analytically, as 
we show in subsection A. For an arbitrary magnetic flux, 
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we discuss a numerical solution in subsection B. 

Consider the BdG equations for spin singlet pairing 



E n v n (r) = 



1 

2 m 



ihV 



— (i 
2m V 



hV 



A(r) - M 



: A(r) 



u n (r)+ Au„(r) 

5 

w n (r)+ A*u n (r) 
(15) 



with the self-consistency condition (gap equation) for the 
order parameter A(r): 



A(r) = K^u n (r)<(r)tanh 



En 

2T 



(16) 



where V is the local pairing potential. For an annulus 
of finite width we separate the angular part of the quasi- 
particle wave functions u n (r), v n (r) using polar coordi- 
nates r = (r, ip) and the ansatz 



u B (r,<p) = u n (r)e^ k+ ^ 
v n {r^) = v n {r)e^ k -^ 



(17) 



where k and q arc either both even or both odd integers. 
Thus hk is the angular momentum as for the II? ring and 
n = (k, p) with the radial quantum number p. The order 
parameter factorizes into A(r,ip) = A(r)e lqip where the 
radial component 



A(r) = T/ ^ u n (r)v n (r) tanh 



2T 



(18) 



is real. For a magnetic flux $ threading the interior of 
the annulus we choose the vector potential A(r, Lp) — 
e v <I>/(27rr), where e v is the azimuthal unit vector. With 
(j) = (e/hc)Q and 



-d r {rd r ) + ^{-id v ±c}>) 2 (19) 




FIG. 5: Annulus with inner radius Ri and outer radius R.2- 
For a magnetic flux threading the interior of the annulus, the 
radial part of the Bogoliubov - de Gennes equations is solved 
numerically with a discretized radial coordinate. 



the BdG equations therefore reduce to radial differential 
equations for u n (r) and v n (r): 



E n u n (r) 



E n v n (r) 



r ffl g fi2/2 

^-(^)-f\ 
Am r Zmr^ 

2;2 



r H 2 d r 
2m r 



(rd r )- 



2mr 2 



-[i u n (r)+ A(r)v n (r) 

1 

v n (r)+ A(r)u n (r) 
(20) 



with the canonical angular momenta 



hi,, 



(k 



til v = -{k-q + 24>). 



(21) 
(22) 



The number q plays the same role as in the previous 
section. Here we choose q for each value of the flux to 
minimize the total energy of the system. The flux for 
which q changes to the next integer can therefore deviate 
from the values (2n — l)/4, where we fixed the change of 
q for the ID model. 



Hankel-Function Ansatz 



A natural choice of an ansatz for the solutions of the 



coupled differential equations ( 20 1 are linear combina- 



tions of the Hankel functions flj and Hj 2 \ since they 



are individually solutions of the uncoupled equations ( 20 1 
for A(r) = 0: 



l dr (rd r ) - ^) ffpV) = 7 2 #, (1 ' 2 V)- 
We therefore take u n (r) and v n (r) of the form: 
u a {r)=U B k W (7» + c^ 2) ( 7 » 



v n {r) 



In 



In 

(2), 



(23) 

(24) 
(25) 



The equations ( 20 1 then become 

h 2 



E n u n (r) 



E n v n (r) 



2m 



MY + m 



u n (r) + A(r)v n (r) 



r s2 



2m 



(26) 



v n (r) + A(r)u n (r) 



The coefficients 7^ and c„ with a — u, v are fixed by 
the open boundary conditions, for which w n (r) and v a {r) 
vanish on the inner and outer boundaries of the annulus: 
Un(-Ri) = u n (R 2 ) = and v a (R{) = v n (R 2 ) = 0. This 
generates the defining equations for 7^ and 



< 1} (7gfii) 
< 2) (7^i) 



< 1) (7g^2) 



(27) 
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For all integer and half-integer values of flux, q = 2<f> 
in the ground state, thus l u = l v = k/2. Assuming a 
constant order parameter A(r) = A, the r-dependcncc 
drops out from Eqs. (26) and we find the eigenvalues and 



eigenvectors of the usual BCS type 



EL 



'h 2 

2m 



+ A2 



(28) 



with 7„ = 7^ = 7^ and 



1 

u„ = - 
2 


1 + 


2 

\2m ln 


bp 


/En 


1 

Vn = - 
" 2 


1 - 


(K 2 2 
^2m 7n 


bp 


/E n 



(29) 
(30) 



These are the two distinct classes of superconducting 
states as discussed for the ID loop: for integer flux val- 
ues, A is given by summing over all even angular mo- 
menta k, whereas for half-integer flux values, A is ob- 
tained by summing over odd angular momenta. 

For general values of magnetic flux, l u and l v are 
different and so are 7^ and 7^. The r-dependence of 
Un(r) is therefore different from u n (r) as contained in 
Eqs. ( 24|25 ). In App. A we analyze the solution of the 
uncoupled Eqs. (20 1 for A = and find that the eigen- 



functions account for the flux induced Doppler shift by 
shifting their nodes closer together or further apart — 
most importantly, u n (r) shifts its nodes in the opposite 
direction than does v n (r). This implies that u n (r) and 
v n (r) with the ansatz of Eqs. (24 1 and (25 1 cannot be 
solutions of the coupled Eqs. (26 1 for A(rj=^ 0. 

Moreover, we show in App. A for the limit of a thin 
annulus (i?i ^> R 2 — Ri) that both the Doppler shift and 
the shift of the nodes of u a (r) and v n (r) are in leading 
order linear functions of q—2(j). It is therefore not possible 



to find an approximate solution of Eqs. ( 20 I that contains 
the effects of the Doppler shift but neglects the shift of 
the nodes. Consequently, we have to resort to a numerical 
solution of the radial component of the BdG equations. 



B. Self-Consistent Numerical Solution 



The numerical solution of Eqs. (|20j) is achieved by 
discretizing the interval [Ri,R 2 ] for the radial coordi- 
nate r into M radii r,, which defines the grid constant 
a = (R2 — R\)/M. In this way we obtain for each 
angular momentum hk M radial eigenstates (channels), 
which correspond to the M eigenstates with the lowest 
eigenenergies E n of the continuum model. On this set 
of M radial coordinates, we use the symmetric discrete 
differential operators d l f(r i ) = [f(r i+1 ) - f(ri-i)]/a 
and d?f(n) = [/(r j+1 ) + /(r«_i) - 2/(r i ]]/a 2 . Insert- 
ing these discrete operators into Eqs. ( |20[ ) and using 
(l/r)d r rd r = (l/r)d r + d 2 , one obtains the eigenvalue 



equation 



t + pt 
A 



A 

-i-pt 



E„ 



(31) 



where u n and v n are real and the operators i, p 1 ^, and A 
are defined through 

iu n (n) = t[u n (r i+ i) + u n (n_i)] 

+ t— [u n (r i+1 ) - u„(ri_i)], (32) 

n 



and 



Au n (r l ) = A(n)u n (ri) 



-r - 2 

2 'a z 



(33) 
(34) 



where t = h 2 /{2ma 2 ). A self-consistent solution of 
Eq. (31) and the gap equation 



A(r;) = Vy^^u n (ri)v n (ri) tanh 



E^ 
2T 



(35) 



is found iteratively. The operator i consists of a sym- 
metric and an antisymmetric part with respect to 7"i_i 
and rjA-i. In order to ensure that the eigenvalues of 
Eq. (31 1 are real, the prefactor of the second, anti- 



symmetric term in Eq. ( 32 ) must be smaller or equal 



to the prefactor of the symmetric term, which means 
M > (R 2 - R 1 )/2R 1 . This condition is fulfilled since 
M = (R 2 - R^/a > (R 2 - Rx)/2Rx. 

Once the eigenfunctions of Eq. (31 1 are known, we ob- 



tain the current by evaluating the expectation value of 
the gauge invariant current operator^. The expectation 
value J(r) of the circulating current is found using a Bo- 
goliubov transformation and the ansatz (17 1 in polar co- 
ordinates: 



J(r) 



he 



E 



mr)f(En) - J^(r)f(-E n )} , (36) 



with 



J n » =^Im 
m 



'-d^ ot a (r,<p) 

r r 



= a n (r) 

m r 



(37) 



for a = u, v. The contribution of each quasi-particle state 
to the total current is therefore determined by its angu- 
lar velocity l a . The radial quantum number p and the 
A-dependence enter only through the occupation proba- 
bility which is controlled by the cigencnergy E n . Further, 
the total energy of the system is given by 



E = if E E " E K(ri)f(E n ) + vl(r t )f(-E n )} 



(38) 
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FIG. 6: Non-self-consistent calculation of current and energy 
at T = 0. The circulating current (upper panel) in an annulus 
with an inner radius i?i = 100a and an outer radius R2 = 
150a is shown for fixed, (^-independent A = (light blue 
line), A = 0.002* (blue line), A = 0.004* (dark blue line), 
A = 0.006* (black line). The lower panel shows the difference 
between the total energy of the annulus as a function of (f> 
and the total energy at zero flux for the same values for A as 
above. 



C. Results 

The results of the non self-consistent calculations for 
the circulating current and the total energy at T = 
and for fixed values of A are displayed in Fig. [6] In the 
normal state (A = 0), there are ~ M eigenstates close 
enough to Ep to cross Ep as a function of <f>, unlike in 
small ID rings where only one state crosses Ep. For 
each crossing, a small jump appears in the current as 
a function of <f>. There is a larger jump at the value 
of <j) where the energies of the even-g and odd-g states 
become degenerate and q switches to the next integer. 
The shape of this function depends on the distribution 
of eigenenergies close to Ep and therefore on microscopic 
details of the geometry of the annulus and the Fermi 
energy Ep. A finite A allows for a flux regime with direct 
energy gap and no crossings of Ep, thus in this regime 
the current is linear and the total energy quadratic in 
<f>. For the largest value (A = 0.006*) shown, there is a 
direct gap for all values of <p. Even for this value of A, 
the current and the energy are not exactly /ic/2e-periodic 
because of the energy diffe rence of the even and odd q 
states in finite systems^^il. 




-1/2 -1/4 1/4 1/2 





FIG. 7: Self-consistent calculations for the same annulus as in 
Fig. [6] In addition, the top panel displays the self-consistent 
order parameter A as a function of cf>. The lines correspond to 
the pairing interaction V — (light blue line), V = 0.28* (blue 
line), V = 0.32* (dark blue line), V = 0.38* (black line). The 
black arrows mark the positions of the g-jump for V = 0.38* 
and V = 0.32*. 



The introduction of self-consistency in A does not fun- 
damentally change these basic observations (Fig. |7| . The 
crossover is then controlled by the pairing interaction 
strength V, for which we chose such values as to repro- 
duce the crossover from the normal state to a state with 
direct energy gap for all flux values. The order parame- 
ter A is now a function of 4>. If A(0 = 0) < 0.006* (cf. 
Fig. [6]), the gap closes with <$> and A decreases whenever 
a state crosses Ep. At these flux values we observe a 
sharp increase in the total energy of the annulus. Unlike 
in ID, A does not drop to zero at the closing of the en- 
ergy gap, but decreases stepwise. In two or three dimen- 
sions, A remains finite beyond <p c because it is stabilized 
by contributions to the condensation energy from pairs 
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with relative momenta perpendicular to the direction of 
the current flow and the closing of the indirect energy 
gap does not destroy superconductivity^ 12 ^. Apart from 
these steps, the current (energy) shows the standard lin- 
ear (quadratic) behavior. 

The offset of the g-jump is only relevant for values of 
V for which A is finite for all <fr. In Fig. [7j the offset 
is clearly visible for the largest two values of V (marked 
with black arrows). Its sign depends on the geometry of 
the annulus and the pairing interaction V — the offset 
changes sign for increasing V (cf. Ref . |2"T|) . 

Experimentally more relevant is to control the 
crossover through temperature. With the pairing inter- 
action V sufficiently strong to produce a T = energy 
gap much larger than the maximum Doppler shift, the 
crossover regime is reached for temperatures slightly be- 
low T c . For the annulus described in Fig. [8] the crossover 
proceeds within approximately one percent of T c . The 
crossover regime gets narrower for larger rings propor- 
tional to the decrease of the Doppler shift. In the limit 
of a quasi ID ring of radius R we can be more precise: If 
we define the crossover temperature T* by A(T*) = A c 
and assuming A c <C A, we can use the Ginzburg-Landau 
form of the order parameter 



ACT) 
A(0) 



1.75a /l 



T 



and obtain 
T n — T* 



A' 2 



El 



3.1A(0) 2 12.4A(0) 2 i? 2 3.1T2R 2 ' 



(39) 



(40) 



For a ring with a radius of 2500 lattice constants (« 10 
/im) and A(0) = O.Olt (w 3 meV) one finds the ratio 
(T c - T*)/T c « 1.3 x 10~ 4 . This is in reasonable quali- 
tative agr eeme nt with the experimental results of Little 
and ParkiPESl, discussed also by TinkhanP3. Their the- 
oretical prediction is similar to Eq. ( 40 1 , up to a factor 



in which they include a finite mean free path. Moreover, 
they do not include the difference introduced through 
even and odd q states. This difference was considered in 
calculations of T c by Bogachek et aZP^in the one-channel 
limit. In Eq. (40) the value of A(0) is in fact different 
for even and odd q. Although quantitative predictions 
of T c — T* of the theory presented here might be too 
large compared to the experiment; it serves as an upper 
limit, because it describes the maximum possible persis- 
tent current. Scattering processes in real systems will 
further reduce T c — T* . 

For temperatures close to T c , the difference of the 
eigenenergies of even and odd q states is less important 
than at T = 0. Thus the deviation from the hc/e- 
periodicity of the current and of the order parameter 
is smaller. Furthermore, persistent currents in the nor- 
mal state are exponentially small compared to the per- 
sistent supercurrents below T c . Their respective hc/e- 
periodic behavior is therefore essentially invisible in the 
flux regime where A — 0. For the annulus described in 



0.016 




FIG. 8: The order parameter A and the persistent current 
for the temperature driven transition from the normal to the 
superconducting state in an annulus with inner radius Ri = 
30a and outer radius R2 = 36a. The pairing intercation is 
V — 0.7 1, with a critical temperature of T c w 0.0523 1 for 
zero flux. For these parameters A(T = 0) O.li. The lines 
correspond to the temperatures T = 0.0513t (black line), T = 
0.0520t (dark blue line), T = 0.0522t (blue line). Notice that 
A is slightly different for the flux values <f) — and = ±1/2. 



Fig. U the difference between A(</> = 0) and A(<f> = 1/2) 
is still visible, but the corresponding differences in the 
current are too small. 



IV. 



CONCLUSIONS 



We have described the crossover from the /ic/e-periodic 
persistent currents as a function of magnetic flux in a 
metallic loop to the ft.c/2e-periodic persistent supercur- 
rent in a ID loop as well as in a multi-channel annulus. 
While a, ID superconducting ring is a rather idealized sys- 
tem, it proves valuable for discussing the physics of this 
crossover. A ring with a radius smaller than half the su- 
perconducting coherence length, shows an /ic/e-periodic 
super current, which reaches the critical current at a crit- 
ical flux value <p c , determined by the flux dependent clos- 
ing of the gap. Assuming that this relation remains un- 
changed on a ring with finite thickness d <C R, as indeed 
suggested by the multi-channel model, R c would be of 
the order of 1/im for aluminum rings. In two or three di- 
mensions, A remains finite beyond <f> c . The temperature 
controlled crossover, while cooling through T c , appears 
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within a temperature window proportional to 1/R 2 and 
thus appears hard to detect in experiment. 
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APPENDIX A: DOPPLER SHIFT AND NODES 
OF THE HANKEL FUNCTION ANSATZ 

The ansatz for u n (r) and v n (r) with two Hankel 
functions (Eqs. (24 1 and ( 25 1 ) solves the normal state 
Schrodinger equation for the annulus as well as the BdG 
equations in the superconducting state with integer and 
half-integer flux values. In this appendix we show that 
it is not possible to construct an approximate analytic 
solution for the superconducting annulus that includes 
the effect of the Doppler shift. In this case u n (r) and 
j)„(r) have the independent eigenenergies (H 2 /2m)^ 2 
and (S72m) 7 f. 

For this purpose we analyze the relation between the 
Doppler shift of the eigenfunctions of the annulus in the 
normal state (A = 0) and the shift of their nodes with re- 
spect to the radial coordinate, using the following asymp- 
totic form for the Hankel functions 28 : 



H, 



(1,2) 



7r I tan x 



exp 



±i ( I tan x 



(Al) 



which approximates Hj'^ 2 ' for I 1. Choosing x — 
acos(l/jr) leads with tan(arccos x) — yl — x 2 /x to 



(1,2) 



H, 



x exp 




-1/4 



2 I 

— — I arccos — 



7r 



(A2) 



(1 2) 

Thus Eq. (A2| approximates H\ ' (77") for jr ^> 1. In- 
serting Eq |A2| into the boundary conditions (27 1 de- 



termines the constants and 7^: 



exp 



2z - J)] = exp [2i (DZ{R 2 ) \ 



with 



D«{r) = r\U 1 Z) 2 -^-l a arccos 



/ 2 



7n r 



(A3) 



(A4) 



The wave functions u n (r) and v n (r) (Eqs. ( 24|25 1) be 



u n (r) = u n \ 



7n r 



-1/4 



e i[Dim+l] S i n [Dl{ r ) - Dl(Rx)\ , (A5) 



v n (r) 



In 



-1/4 



^lsinK(r)-TO)]. (A6) 



The vanishing of the wavefunction for r = R2 therefore 
implies that 



D%(R 2 ) - D^(Ri) = -TTp 



(A7) 



for an integer p, which determines 7^. In the limit of a 
thin annulus (i?i 3> R2 — Ri), we expand D^(r) in 1/r 
and find 



D^(r) - D^R,) ^ (r - R,) 



p 



7; 



2 7 «ri?! 



(A8) 



With this asymptotic form the boundary condition ( A7 1 
becomes a quadratic equation in 7° : 



a\2 



up 



"7„ 



R\ — R2 2R1R2 
which has the positive solution 



0. 



(A9) 



np 



Ri — R 2 




Ri — R2 



P 



2RiR 2 



(A10) 



This is the simplest possible approximation for the 
eigenenergies of the uncoupled equations (A = 0) of the 
annulus containing the Doppler shift, which is controlled 
by Z 2 . The flux (f> enters l u a,nd l v with different signs 
(see Eq. 22)). Thus, if 7^ decreases as a function of (f>, 
7„ increases. Since q — 2<j) < 1 in the ground state, the 
Doppler shift (h 2 /2m)[^ 2 (q - 2(f) - 7„ 2 (0)] is linear in 
leading order in {q — 2(f)j\fR\. 

We further find the nodes r nm of u n (r) and v n (r) by 
setting expression ( |A8[ ) equal to nm, where m is a posi- 
tive integer, and solving it for r > 0: 



Ri 



P 



2 1 £ 2 R 1 




I 2 nn 
2lZ 2 Ri ~ 7^ 



2/ 2 

r^a2 
In 



(AH) 



The shift of the nodes r nm (q — 2<fi) — r nm (0) as a function 
of flux is again linear in (q — 2<p) / y/R~i to leading order. 



10 



Thus both the Doppler shift and the nodes of u n (r) shift 
linearly with <p and conversely when compared with the 
Doppler shift and the nodes of v n (r). 

The coupled Eqs. (26) for A ^ resulting from the 
ansatz ( |24|25| with non integer (or non half-integer) flux 
can be solved only by wave functions u n (r) and v a (r) 



with the same r-dependence. To obtain a solution of this 
problem, one can expand the wave functions as a sum 
of Hankel functions and numerically solve for the coeffi- 
cients or directly solve the coupled differential equations 
numerically. 
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